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Abstract 

In the framework of the Lindblad theory for open quantum systems we deter- 
mine the degree of quantum decoherence and classical correlations of a harmonic 
oscillator interacting with a thermal bath. The transition from quantum to clas- 
sical behaviour of the considered system is analyzed and it is shown that the 
classicality takes place during a finite interval of time. We calculate also the deco- 
herence time and show that it has the same scale as the time after which statistical 
fluctuations become comparable with quantum fluctuations. 
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1 Introduction 

The transition from quantum to classical physics and classicality of quantum systems 
continue to be among the most interesting problems in many fields of physics, for both 
conceptual and experimental reasons [HOE]. Two conditions are essential for the 
classicality of a quantum system [HE]: a) quantum decoherence (QD), that means 
the irreversible, uncontrollable and persistent formation of a quantum correlation (en- 
tanglement) of the system with its environment [6], expressed by the damping of the 
coherences present in the quantum state of the system, when the off-diagonal elements 
of the density matrix of the system decay below a certain level, so that this density 
matrix becomes approximately diagonal and b) classical correlations (CC), expressed 
by the fact that the Wigner function of the quantum system has a peak which follows 
the classical equations of motion in phase space with a good degree of approximation, 
that is the quantum state becomes peaked along a classical trajectory. The necessity 
and sufficiency of both QD and CC as conditions of classicality are still a subject of 



debate. Both these conditions do not have an universal character, so that they are not 
necessary for all physical models. An important role in this discussion plays the temper- 
ature of the environment and therefore it is worth to take into account the differences 
between low and high temperature regimes. For example, purely classical systems at 
very high temperatures are described by a classical Fokker-Planck equation which does 
not follow any trajectory in phase space (for very small kinetic energy, compared to the 
thermal energy, when the probability distribution becomes essentially independent of 
momentum), so that in this case CC are not necessary. Likewise, one can have a classi- 
cal behaviour if the coherences are negligible, without having strong CC (for example, 
in the case of a classical gas at finite temperature) and the lack of strong correlations 
between the coordinate and its canonical momentum does not necessarily mean that 
the system is quantum. On the other hand, the condition of CC is not sufficient for a 
system to become classical - although the Wigner function can show a sharp correlation 
in phase space, the quantum coherence never vanishes for a closed system which has a 
unitary evolution. Likewise, in the low temperature quantum regime one can observe 
strong CC. For example, in the case of a purely damped quantum harmonic oscillator 
(at zero temperature), the initial coherent states remain coherent and perfectly follow 
classical trajectories of a damped oscillator, but CC are not sufficient for classicality. 

In the last two decades it has became more and more clear that the classicality 
is an emergent property of open quantum systems, since both main features of this 
process - QD and CC - strongly depend on the interaction between the system and 
its external environment [31 d [HI [91 [10l HH [12l [131 LI31 LTS] . A remarkable aspect of 
the current research helping in understanding the nature of the quantum to classical 
transition is that for the first time there have recently been carried on experiments 
probing the boundary between the quantum and the classical domains in a controlled 

way pn El mi eh na HQ. 

The role of QD became relevant in many interesting physical problems of field the- 
ory, atomic physics, quantum optics, quantum information processing, quantum gravity 
and cosmology, and condensed matter physics. We mention here only a few of these 
problems [21 [3] : to understand the way in which QD favorizes the quantum to classical 
transition of density fluctuations; to study systems of trapped and cold atoms (or ions) 
which may offer the possibility of engineering the environment, like trapped atoms inside 
cavities, relation between decoherence and other cavity QED effects (such as Casimir 
effect); on mesoscopic scale, decoherence in the context of Bose-Einstein condensation. 

In many cases one is interested in understanding the specific causes of QD just 



2 



because one wants to prevent decoherence from damaging quantum states and to protect 
the information stored in quantum states from the degrading effect of the interaction 
with the environment. Thus, decoherence is responsible for washing out the quantum 
interference effects which are desirable to be seen as signals in some experiments. QD 
has a negative influence on many areas relying upon quantum coherence effects, such 
as quantum computation and quantum control of atomic and molecular processes. In 
the physics of information and computation, decoherence is an obvious major problem 
in the implementation of information-processing hardware that takes advantage of the 
superposition principle [22\. 

In most of literature, QD has been studied for a system coupled to an environment 
or thermal bath with many degrees of freedom. The main purpose of this paper is to 
study QD and CC for a harmonic oscillator interacting with an environment in the 
framework of the Lindblad theory for open quantum systems. More concretely we 
determine the degree of QD and CC and the possibility of simultaneous realization 
of QD and CC for a system consisting of a harmonic oscillator in a thermal bath. 
For that purpose, we first find the evolution of the density matrix and of the Wigner 
function of the considered system and then we apply the criterion of QD and CC. We 
consider different regimes of the temperature of environment. It is found that the system 
manifests a QD which increases with time and temperature, whereas CC are less and 
less strong with increasing time and temperature. 

The organizing of the paper is as follows. In Sec. 2 we review the Lindblad 
master equation for the damped harmonic oscillator and in Sec. 3 we derive the master 
equation in coordinate representation and the corresponding Fokker-Planck equation in 
the Wigner representation and determine the density matrix and Wigner function of 
the considered system. Then in Sec. 4 we investigate QD and CC and analyze them 
quantitatively. In Sec. 5 we calculate the decoherence time of the system and discuss 
the transition from quantum mechanics to classical statistical mechanics. A summary 
and concluding remarks are given in Sec. 6. 

2 Lindblad master equation for the harmonic oscil- 
lator 

Here we review the Lindblad's axiomatic formalism based on quantum dynamical semi- 
groups. The irreversible time evolution of an open system is described by the following 
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general quantum Markovian master equation for the density operator p(t) |23 | [241125] : 
^ = ~[H lP (t)] + k^ i[Vjp{t) ^ ] + M>P®V}])- C 1 ) 

H is the Hamiltonian of the system and Vj, Vj are operators on the Hilbert space of H, 
which model the environment. In order to obtain, for the damped quantum harmonic 
oscillator, equations of motion as close as possible to the classical ones, the two possible 
operators V\ and V2 are taken as linear polynomials in coordinate q and momentum p 
[2Hj |2"TI |2"8"] and the harmonic oscillator Hamiltonian H is chosen of the general quadratic 
form 

£*/ x 1 9 moo 2 9 

H = H + ^(qp + pq), H = —p 2 + —q 2 . (2) 

With these choices the master equation ([Tj) takes the following form [271 EE]: 

-77 = —AH Q ,p\ - — (A + p)[q,pp + pp] + ^(A-p)[p,pg + gp] 

-3»[ g> [g,p]] - [p,p]] + ^([g, [p,p]] + [p, [g,p]]). (3) 

The quantum diffusion coefficients D pp , D qq , D pq and the dissipation constant A satisfy 
the following fundamental constraints [27J EE]: D pp > 0, D qq > and 

D PP D qq - D 2 pq > ^f. (4) 

In the particular case when the asymptotic state is a Gibbs state pc(oo) = e~kr /Tre - ^ , 
these coefficients become |27l [28] 



A + p too X — p Ti too 

Dr>v= amuj coth — — -, D aa = coth — — , D va = 0, (5) 

pp 2 2fcT' 99 2 mcj 2kT' pq K J 

where T is the temperature of the thermal bath. In this case, the fundamental con- 
straints are satisfied only if A > p and 

(A 2 - P 2 ) coth 2 > A 2 . (6) 

From the master equation ([3]) we obtain the following equations of motion for the 
expectation values of coordinate and momentum J27J [28] : 



-(A-p)^(t) + -a p (t), (7) 
m 

-mu 2 a q {t) - (A + p)a p {t). (8) 



dt q m 

dcr p (t) 
dt 
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In the underdamped case {u > //) considered in this paper, with the notation Q 2 = 
ui 2 - /j, 2 , we obtain [271 EH]: 

a q (t) = e~ A *((cosfit + ^sinttt)a q (0) + sin Qta p {0)), (9) 

i U II ii U 

9 

TTl/UJ Li 

(cos fit — — sinfit)crp(O)) (10) 

and (Tg(oo) = 0p(oo) = 0. 

Lindblad has proven [26J that in the Markovian regime the harmonic oscillator mas- 
ter equation which satisfies the complete positivity condition cannot satisfy simultane- 
ously the translational invariance and the detailed balance (which assures an asymptotic 
approach to the canonical thermal equilibrium state). The necessary and sufficient con- 
dition for translational invariance is A = /z [26l [271 EH] • In this case the equations of 
motion (JTj) and jSJ) are exactly the same as the classical ones. If A ^ //, then we violate 
translational invariance, but we keep the canonical equilibrium state. 

The relation (HI) is a necessary condition for the generalized uncertainty inequality 

h 2 

<r qq (t)a PP (t) - a 2 pq (t) > — (11) 

to be fulfilled, where a qq and a pp denote the dispersion (variance) of the coordinate and 
momentum, respectively, and a pq denotes the correlation (covariance) of the coordinate 
and momentum. The equality in relation (fTTl) is realized for a special class of pure 
states, called correlated coherent states [29] or squeezed coherent states. 

The asymptotic values a qq (oo), a pp (oo), a pq (oo) do not depend on the initial values 
a qq (Q), <t pp (0), cr pq (0) and in the case of a thermal bath with coefficients (jSJ), they reduce 
to [27JEH] 

h hu Timu Too 

aqq ^ = 2^ 2fcT' a Pp(°°) = ~Y~ 2kT' o Pq {°°) = 0- (12) 

In the following, we consider a general temperature T, but we should stress that 

the Lindblad theory is obtained in the Markov approximation, which holds for high 

temperatures of the environment. At the same time, the semigroup dynamics of the 

density operator which must hold for a quantum Markovian process is valid only for the 

weak-coupling regime, with the damping A obeying the inequality A <C uj. 

3 Density matrix and Wigner distribution function 

We consider a harmonic oscillator with an initial Gaussian wave function 
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where cr qq (0) is the initial spread, u pq (0) the initial covariance, and cr g (0) and cx p (0) 
are the initial averaged position and momentum of the wave packet. The initial state 
( IT3l) represents a correlated coherent state [29] with the variances and covariance of 
coordinate and momentum 

hS / \ hmu) , „ hr , 

'-(») = SSJ. MO) = ^372). ^(°» = VT=7I' (14) 

Here, 5 is the squeezing parameter which measures the spread in the initial Gaussian 
packet and r, with |r| < 1 is the correlation coefficient at time t — 0. The initial 
values (fl4"|) correspond to a minimum uncertainty state, since they fulfil the generalized 
uncertainty relation 

- <&(o) = x- (15) 

For 5 = 1 and r = the correlated coherent state becomes a Glauber coherent state. 
For a given temperature T of the bath and for any parameters 5 and r the inequality 
(JIJJ) alone determines the range of values of the parameters A and fi [30J . 

From Eq. Q we derive the evolution equation in coordinate representation: 

dp ih d 2 d 2 imuj 2 2 /2 
m = ^d^ 2 ~dq ;2 ^^ ~ 
-i(A + - A), + i(A - + MjL + A) + 2]p 

-^(9 - 9')V + B M (| + - 2iD„R(, - + |> (16) 

and in Refs. [211 E21 E3] we transformed the master equation for the density operator 
into the following Fokker-Planck-type equation satisfied by the Wigner distribution 
function W(q,p, t) : 

dW pdW 2 dW ^d, tt7-\ /> J,„. 

1 muo 2 q— — + (A + + (A - jt) — (gWf) 



+ + 2^pZHr. (17) 



m <9g dp dp dq 

^-K + D ^^ + 2D d ' — 
<9p 2 99 <9g 2 pq dpdq 

The first two terms on the right-hand side of both these equations generate a purely 
unitary evolution. They give the usual Liouvillian evolution. The third and forth terms 
are the dissipative terms and have a damping effect (exchange of energy with envi- 
ronment). The last three are noise (diffusive) terms and produce fluctuation effects in 
the evolution of the system. D pp promotes diffusion in momentum and generates deco- 
herence in coordinate q: it reduces the off-diagonal terms, responsible for correlations 
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between spatially separated pieces of the wave packet. Similarly D qq promotes diffusion 
in coordinate and generates decoherence in momentum p. The D pq term is the so-called 
"anomalous diffusion" term. It promotes diffusion in the variable qp + pq, just like both 
the other diffusion terms, but it does not generate decoherence. 

In the high temperature limit, quantum Fokker-Planck equation ffTTj) with coeffi- 
cients ((5]) becomes classical Kramers equation (D pp — > 2m\kT for A = /i) [33] . 

The density matrix solution of Eq. (fT6l) has the general form of Gaussian density 
matrices 

-w^^ - + - - + - < 18 > 

where 

<r(t) = <r M (*)M*) - a 2 pq (t) (19) 
is the determinant of the dispersion (correlation) matrix 

CTqqit) <7 pq {t) 

it) a pp {t) 



(20) 



and represents also the Schrodinger generalized uncertainty function [30] . 

For an initial Gaussian Wigner function (corresponding to a correlated coherent 
state (TT31) ) the solution of Eq. (|T7|) is 

1 i 

W{q,p,t) = ==exp{-— — Aa w {t){q- a q {t)f + a qq {t)(jp - a p (t)Y 

2nJa{t) 2o W 

-2a pq (t)(q-a q (t))(p-a p (t))}}. (21) 



In the case of a thermal bath we obtain the following steady state solution for 

2fcT' 



oo (we denote e = 



< >= <^h^ ex ^-4^[^ + ( * - coth * (22) 

In the long time limit we have also 

W 00 (q,p)= 1 exp{ - — 1— H 2 + -^]}. (23) 
7m coth e n coth e mcu 

Stationary solutions to the evolution equations obtained in the long time limit are 
possible as a result of a balance between the wave packet spreading induced by the 
Hamiltonian and the localizing effect of the Lindblad operators. 
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4 Quantum decoherence and classical correlations 



As we stated in the Introduction, one considers that two conditions have to be satisfied in 
order that a system could be considered as classical. The first condition requires that the 
system should be in one of relatively permanent states (states that are least affected by 
the interaction of the system with the environment, called by Zurek "preferred states" in 
the environment induced superselection description [21 E]) and the interference between 
different states should be negligible. This implies the destruction of off-diagonal elements 
representing coherences between quantum states in the density matrix, which is the QD 
phenomenon. An isolated system has an unitary evolution and the coherence of the 
state is not lost - pure states evolve in time only to pure states. The loss of coherence 
can be achieved by introducing an interaction between the system and environment: an 
initial pure state with a density matrix which contains nonzero off-diagonal terms can 
non-unitarily evolve into a final mixed state with a diagonal density matrix during the 
interaction with the environment, like in classical statistical mechanics. 

The second condition requires that the system should have, with a good approxima- 
tion, an evolution according to classical laws. This implies that the Wigner distribution 
function has a peak along a classical trajectory, that means there exist CC between 
the canonical variables of coordinate and momentum. Of course, the correlation be- 
tween the canonical variables, necessary to obtain a classical limit, should not violate 
Heisenberg uncertainty principle, i.e. the position and momentum should take reason- 
ably sharp values, to a degree in concordance with the uncertainty principle. This is 
possible, because the density matrix does not diagonalize exactly in position, but with 
a non-zero width, i.e. it is strongly peaked about q — q' and very small for q far from 



Using new variables S = (q + q')/2 and A = q — q', the density matrix (fl8l) can be 
rewritten as 



p(E,A,t) 




with the abbreviations 



1 



a{t) 



ha qq {t) 



(25) 



a = 



2a qq (t) 



7 = 



2h 2 a qq {t) 



and the Wigner transform of the density matrix ( J24l) is 



W(q,p,t) 



2nh 



1 




exp{ 



[h/3(q - a q (t)) - (p - cr p (t))] 2 



a(q-a q (t)) 2 }. (26) 



S 



a) Degree of quantum decoherence ( QD 

The representation-independent measure of the degree of QD [¥J is given by the 
ratio of the dispersion l/y/2j of the off-diagonal element p(0, A, t) to the dispersion 



(27) 



2/ a of the diagonal element p(£,0,£) : 



'QD 




which in our case gives 



S QD (t) = —==. {2i 
2Mt) 



The finite temperature Schrodinger generalized uncertainty function (I19p . calcu- 
lated in Ref. [301, nas the expression 

"(t) = ^{e- 4A, [l - (6 + S{1 l _ r2) ) cothe + coth 2 £ ] 

^ - w^-)^ + X 1 ^" ' ' + coth2 E} (29) 

In the limit of long times Eq. (1291) yields 

h 2 



r 



so that we obtain 



which for high T becomes 



o-(oo) = ^- coth 2 e, (30) 



^qd(oo) = tanh^, (31) 



<W°°) = (32) 



We see that Sqd decreases, and therefore QD increases, with temperature, i.e. the 
density matrix becomes more and more diagonal at higher T and the contributions 
of the off-diagonal elements get smaller and smaller. At the same time the degree of 
purity decreases and the degree of mixedness increases with T. 5qd < 1 for T ^ 0, 
while for T = the asymptotic (final) state is pure and 5qd reaches its initial maximum 
value 1. A pure state undergoing unitary evolution is highly coherent: it does not lose 
its coherence, i.e. off-diagonal coherences never vanish. 5qd = when the quantum 
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coherence is completely lost. So, when 8qd = 1 there is no QD and only if 5qd < 1, 
there is a significant degree of QD, when the magnitude of the elements of the density 
matrix in the position basis are peaked preferentially along the diagonal q = q' '. When 
5qd <C 1, we have a strong QD. 

b) Degree of classical correlations ( CC) 

In defining the degree of CC, the form of the Wigner function is essential, but not 
its position around a q (t) and a p (t). Consequently, for simplicity we consider zero values 
for the initial expectations values of the coordinate and momentum and the expression 
( |26l) of the Wigner function becomes 



2 



w ^=^^ JJ %r- afl - (33) 

A ridge of the Wigner function (I33p in phase space is at p = H/3q, showing the correlation 
between q and p. As a measure of the degree of CC we take the relative sharpness of 
this peak in the phase space determined from the dispersion hy/2/y in p in Eq. (|33|) and 
the magnitude of the average of p (p = K/3q) [I]: 

2^/cry 



>cc 



(34) 



where we identified q as the dispersion l/v2a of q. 5cc is a good measure of the 
"squeezing" of the Wigner function in phase space [3]: in the state (TSUI) , more "squeezed" 
is the Wigner function, more strongly established are CC. In the coordinates h/3q — p 
and %(3q (these quantities have the same dimension), 2fi^ and %\(3\/y/a are the lengths 
of the shorter and longer semi-axes of the la contour in phase space and their ratio gives 
5cc- Similarly, in coordinates %(3q — p and q, 2fi v / 7 and 1/yfot are the lengths of the 
shorter and longer semi-axes of the la contour and their product gives the area of the 
la ellipse. We see from Eq. ( 1271) that Sqd is inversely proportional to this area. Besides 
this geometric interpretation, 5qd is also connected with the linear entropy 
For our case, we obtain 



5cc(t) = f^J- (35) 

where a{t) is given by Eq. fl29l) and a pq {t) can be calculated using formulas given in 
Refs. [23 [28]: 

h 1 2uj 2 r 
a pS) = 77^e- 2A '{[M2cothe -5- - n _ ] cos(2fit) 



4ft 2 v 5(1 -r 2 ) ' y/l 



r 



1 1 2u 2 r 
+ujQ(5 - —) sin(2ftt) + fjw(6 + — - 2 coth e) + ^ =}. (36) 



S(l-r*y K ' ^ v 5(1 -r 2 ) - — ^1 



r 
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CTp(t) 




Figure 1: Trajectory in phase space given by the expectation values of coordinate a q (t) 
and momentum cr p (t) for time t G [0, 14] with initial coordinate a q (0) = 6 and mo- 
mentum <Jp(0) = 4 and la contours of the Wigner function, corresponding to an initial 
coherent state (5 = 1) and a squeezed state (5 = 4), for A = 0.2 and ji = 0.1. In all 
figures we use the system of units m = u = h = 1. 



When 5qc is of order of unity, we have a significant degree of classical correlations. The 
condition of strong CC is 5cc ^ 1, which assures a very sharp peak in phase space. 
Since a pq (oo) = 0, in the case of an asymptotic Gibbs state, we get dcci^o) —> °o, so 
that our expression shows no CC at t — > oo. 

c) Discussion with Gaussian density matrix and Wigner function 
We have seen that if the initial wave function is Gaussian, then the density matrix 
( I18p and the Wigner function (121]) remain Gaussian for all times (with time-dependent 
parameters which determine their amplitude and spread) and centered along the trajec- 
tory given by Eqs. ([HD and (fT0~|) . which are the solutions a q (t) and a p (t) of the dissipative 
equations of motion ([7j) and (jSJ). This trajectory is exactly classical for A = \x and only 
approximately classical for not large A — \x. In Fig. 1 there are represented the trajec- 
tory in phase space and two examples of the la contour of the initial Wigner function, 
corresponding to an initial coherent state (5 = 1) and a squeezed state (5 = 4). In 
general, the la contour is defined by the ellipse 

- ^(t)) 2 + M*)(P - o*»(*)) 2 - 2o wW(<? - ^(*))(P " vpit))] = 1- (37) 

In Fig. 1 the center of the ellipses is the point given by the initial expectation values 
(Tq(O) and a q (0). 

To illustrate the dependence on the temperature and time of the degree of QD 
and the degree of CC, we represent them in Fig. 2. We can see that the degree of 
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QD has an evolution which shows that in general QD increases with time and tem- 
perature. We can also see that the degree of CC has a more complicated evolution, 
but the general tendency is that CC are less and less strong with increasing time and 
temperature. 5qd < 1 for non-zero temperature and See is of the order of unity for a 
long enough interval of time, so that we can say that the considered system interacting 
with the thermal bath manifests both QD and CC and a true quantum to classical 
transition takes place. Dissipation promotes quantum coherences, whereas fluctuation 
(diffusion) reduces coherences and promotes QD. The balance of dissipation and fluc- 
tuation determines the final equilibrium value of 8qd- The quantum system starts as a 
pure state, with a Wigner function well localized in phase space (Gaussian form). This 
state evolves approximately following the classical trajectory (Liouville flow) in phase 
space and becomes a quantum mixed state during the irreversible process of QD. 

The squeezing and the correlation of the initial state play also a role in the degree 
of QD and CC. When the squeezing parameter 5 is increasing, both Sqd and Sec are 
decreasing, therefore the squeezing favorizes both QD and CC. Likewise, the increasing 
of the correlation coefficient r leads to the decreasing of the degree of QD, so that when 
we have a larger initial correlation, then QD is stronger. At the same time, the variation 
of the correlation coefficient has a small influence on the general pulsatory behaviour 
of the degree of CC. We also remark from Eq. (13T1) that the asymptotic value of the 
degree of QD does not depend on the initial squeezing and correlation, it depends on 
temperature only. 

In Figs. 3 and 4 we represent the density matrix in coordinate representation ffl8l) 
and the Wigner function ( 12TI) at the initial and final moments of time. The values of 
the density matrix along the diagonal q = q' represent the probability of finding the 
system in this position, while the off-diagonal values represent the correlations in the 
density matrix between the points q and q'. The asymptotic Wigner distribution has an 
axial symmetry, reflecting quantum equipartition. For simplicity, in Figs. 3 and 4 we 
consider zero values for the initial expectations values of coordinate and momentum, so 
that both the density matrix and Wigner function are centered in origin. Of course, as 
we stated earlier, for non-zero initial expectations values of coordinate and momentum, 
the density matrix (Ti~8l) and the Wigner function ( T2T1) are centered along the trajectory 
given by Eqs. ([9]) and (jTUjl . like in Fig. 1. 

From expressions ff27|) and fl34l) we notice that the key parameter which describes 
QD and CC is 7. This coefficient determines the spread of the Wigner function ( 1261) 
around the path in phase space and measures the contribution of non-diagonal terms 
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Figure 2: a: Degree of quantum decoherence 5qd and b: degree of classical correlations 
Sec as functions on temperature T (through C = coth Jj^) and time t for A = 0.2, /z = 
0.1,<J = 4,r = 0. 
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Figure 3: Density matrix p in coordinate representation for A = 0.2, p, = 0.1, 5 = 4, r = 
0; a: |p| at the initial time t = 0; b: at time i — > oo for C = 3; c: for C = 20. 
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in the density matrix (1241) . Therefore, when decoherence increases, the correlations 
between the canonical variables of coordinate and momentum decrease. The extreme 
limit of QD (7 — > 00) is incompatible with CC and that of CC (7 — > 0) is incompatible 
with QD. Their simultaneous realization is not a trivial task: QD requires interaction 
with an environment, which inevitably suppresses CC and produces fluctuations in the 
evolution of the system, whereas classical predictability requires these fluctuations to 
be small. Therefore the existence of the environment is crucial for the quantum to 
classical transition and, consequently, classicality is an emergent property of an open 
quantum system. The shown figures confirm the presence of the relative competition 
which appears between QD and existence of CC, since decoherence (diagonalizing or 
the decreasing of the width of the density matrix) implies a spreading of the Wigner 
distribution function (which is the Fourier transform of the density matrix) along the 
trajectory in phase space, whereas CC require the existence of sharp peaks in the Wigner 
function. Although there exists this competition, there is a broad compromise regime in 
which QD and CC can hold well simultaneously. If the density matrix is not diagonal, 
but the Wigner function becomes peaked along the classical trajectory for long times, we 
do not have, strictly speaking, a classical limit, but only a classical behaviour and CC. 
We regard classical behaviour as a quantum behaviour in which there exists correlations 
between coordinate and momentum or when coherences are negligible, even without 
having CC. 

We can assert that in the considered case classicality is a temporary phenomenon, 
which takes place only at some stages of the dynamical evolution, during a definite inter- 
val of time. Due to the dissipative nature of evolution, the approximately deterministic 
evolution is no more valid for very large times, when the localization of the system is 
affected by the spreading of the wave packet and of the Wigner distribution function. 

In the case of a closed harmonic oscillator (zero-damping limit, A = 0), 5qd = 1 
and the QD phenomenon does not take place. At the same time we obtain (if r = 0) 



Consequently, for an initial coherent state (5 = 1), we see that 5cc — 00 f° r an Y time, so 
that there are no CC. If the initial state is squeezed (5^1), then the system manifests 
CC and it has a classical behaviour. 
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5 Decoherence time scale and transition from quan- 
tum mechanics to classical statistical mechanics 



Diffusion in momentum, which generates the decoherence in coordinate q, occurs at the 
rate set by D pp . In the macroscopic limit, when H is small compared to other quantities 



with dimensions of action, such as \J D pp < (q — q') 2 >, the term in Eq. (IT5|) containing 
Dpp/h 2 dominates and induces the following evolution of the density matrix: 

£ = -£<»-rtV (39) 

Thus the density matrix loses off-diagonal terms in position representation: 

< q\p(t)W >=< q\p(0)W > exp[-^r(q - </)% (40) 

while the diagonal (q = q') ones remain untouched. Quantum coherences decay expo- 
nentially at a rate given by 

D 



f{q~q)\ (41) 



so that the decoherence time scale is of the order of 

h 2 



( 42 ) 



Dpp{q - q 

In the case of a thermal bath, we obtain (see Eq. (J5j)) 

t - 2H (43) 

eco (A + /j)mucr qq (0) cothe' 

where we have taken (q — q') 2 of the order of the initial dispersion in coordinate a qq (0). 

In order to obtain a more precise expression of the decoherence time, we consider 
the coefficient 7 ( J25l) . which measures the contribution of non-diagonal terms in the 
density matrix (pMj) . For short times (At <C 1, fit <C 1), we have: 

2 2 

TflLO V T LOT 

From here we obtain that quantum coherences in the density matrix decay exponentially 
at a rate given by 

2 2 

T r LOT 

2 ^ 6 + j(i- r 2) ) coth e + ~ ^ ( i_ r2 ) ) coth e " A " * ~ 5^T^ ] (45) 
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and then the decoherence time scale is 

tdeC ° = 2[A(5+^)cothe + M5-^)coth e -A-/i-^]- (46) 

The decoherence time depends on the temperature T and the coupling A (dissipation 
coefficient) between the system and environment (through the diffusion coefficient D pp ), 
on the squeezing parameter S that measures the spread in the initial Gaussian packet and 
on the initial correlation coefficient r. We notice that the decoherence time is decreasing 
with increasing dissipation, temperature and squeezing. 
For r = we obtain: 

tdeco = 2(A + //)(* cothe -1) (47) 
and at temperature T = (then we have to take /i = 0), this becomes 

tdeco = 2A(5-1)' (48) 

We see that when the initial state is the usual coherent state (5 = 1), then the decoher- 
ence time tends to infinity. This corresponds to the fact that for T = and 5 = 1 the 
coefficient 7 is constant in time, so that the decoherence process does not occur in this 
case. 

At high temperature, introducing the notation 

2kT 1 (49) 



expression fH6|) becomes 



1 

tdec ° " 2[A(5 + ^) + M5-^)K (50) 

If, in addition r = 0, then we obtain 

%uj , . 

tdeco = 4(X + fi)SkT { } 

In Ref. [30] we studied the behaviour of the generalized uncertainty function a(t) 
( |29l) . For short times we obtained 

ait) = ^{1 + 2[\{8 + jt^w ) cothe + fi(6 - ^ ) cothe - 2\}t}. (52) 

This expression shows explicitly the contribution for small time of: (i) uncertainty that is 
intrinsic to quantum mechanics, expressed through the Heisenberg uncertainty principle 
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and (ii) uncertainty due to the coupling to the thermal environment, which has two 
components, dissipation and diffusion (this last component is responsible for the process 
of decoherence) . From Eq. (l52|) we can determine the time td when statistical (thermal) 
fluctuations become comparable with quantum fluctuations. At high temperature we 
obtain 

1 

U ~ 2r[\(6 + ^) + ^ ~ s^W (53) 

By statistical (thermal) fluctuations we mean the fluctuations produced by diffusion, 
that arise in the generalized uncertainty function a(t) from the coupling of the har- 
monic oscillator to the thermal bath at arbitrary temperature T, even at T = (when 
the diffusion coefficient still has a non-zero value). By quantum fluctuations we mean 
fluctuations of the quantum harmonic oscillator at zero coupling with the thermal bath. 

As expected, we can see that the decoherence time td eco has the same scale as 
the time td after which statistical fluctuations become comparable with quantum fluc- 
tuations. The values of tdeco and td become closer with increasing temperature and 
squeezing. 

When t 3> t re i, where t re i « A^ 1 is the relaxation time, which governs the rate of 
energy dissipation, the particle reaches equilibrium with the environment. Indeed, the 
uncertainty function a(t) ( 1291) is insensitive to A, fi, 5 and r and approaches 

a BE = ^ CQth 2 e _ (54) 

4 

This is the Bose-Einstein relation for a system of bosons in equilibrium at tempera- 
ture T, obtained also in quantum Brownian models for the weak coupling at arbitrary 
temperature. In the case of T = we approach the limit of pure quantum fluctuations, 

h 2 , , 

= -j, (55) 

which is the quantum Heisenberg relation. At high temperatures T (T ^> hu/k) we 
obtain the limit of pure thermal fluctuations, 

a MB = (KTy (56) 

UJ 

which is a Maxwell-Boltzmann distribution for a system approaching a classical limit. 

For all macroscopic bodies the dissipation term becomes important much later after 
the decoherence term has already dominated and diminished the off-diagonal terms. If 
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we compare the time scales of these two terms (decoherence rate (HT|) and relaxation 
rate) we get for high temperatures (for fi = 0) 

Decoherence rate D pp(f q ^ muJ , %uj mkT , ,. 9 

-— = — \ = — (o - q') coth— — « — o-(q-q ) - (57) 

Relaxation rate A 2Ti yH H ' 2kT h 2 y ' y ' 

In most typical situations this is a huge number. For example, for a mass of 1 g at room 
temperature (T = 300 K) and for a separation of q — q 1 — 1 cm, the decoherence time 
scale tdeco is approximately 10 40 times shorter than the relaxation time t re i, so that in the 
macroscopic domain QD occurs very much faster than relaxation. We remark also that 
tdeco can be of the order of t re \ for sufficiently low temperatures and small wave packet 
spread (small squeezing coefficient). At the same time we have to remind that Lindblad 
theory is obtained in the Markovian approximation, when the characteristic time scales 
of the considered processes are larger than the characteristic time scales of the thermal 
bath. When the decoherence time scale is much shorter than the relaxation time and at 
the same time the decoherence becomes faster even than the environment time scales, 
then Markovian approximation is no more valid and the use of non-Markovian quantum 
diffusion models becomes preferable. 

We have seen that a necessary condition for a system to behave classically is the QD 
process. The time scale after which the system has a classical behaviour, with an evolu- 
tion described by a classical probability distribution, is determined by the decoherence 
time tdeco, when the density matrix becomes approximately diagonal very rapidly during 
the system-environment interaction. On the other hand, one often regards the regime 
where statistical fluctuations begin to surpass quantum fluctuations as the transition 
point from quantum to classical statistical mechanics and identifies the high temper- 
ature regime of a system as the classical regime. Above it was shown that these two 
criteria of classicality are equivalent: the time when the quantum system decoheres is 
comparable with the time when statistical fluctuations overtake quantum fluctuations 
(tdeco ~ td). This result is a new confirmation of the previous similar results [36, 37, 38J. 
However the regime after statistical fluctuations dominate should not be called classical. 
After the decoherence time, although the system is describable in terms of probabilities, 
it can not yet be regarded as classical because of the spin-statistics effects. It has to be 
described by non-equilibrium quantum statistical mechanics. After the relaxation time 
the system is correctly treated by the equilibrium quantum statistical mechanics, and 
only at a sufficiently high temperature, when the spin (Fermi-Dirac or Bose-Einstein) 
statistics can be represented by the Maxwell-Boltzmann distribution function, it can be 
considered in a classical regime [301 E3 EE] (see Eqs. (1541) . ( 156]) ). 
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6 Summary and concluding remarks 



We would like to remark that there are no precise quantitative criteria in literature 
for classicality and at the same time there exist some ambiguity and even unclarities 
concerning the characteristics of the quantum to classical transition. In the present 
paper we have studied QD and CC with the Markovian equation of Lindblad in order 
to understand the transition from quantum to classical mechanics for a system con- 
sisting of an one-dimensional harmonic oscillator in interaction with a thermal bath in 
the framework of the theory of open quantum systems based on quantum dynamical 
semigroups. Our results may be summarized as follows. 

(1) Using the criterion of QD for the considered model, we have shown that QD 
in general increases with time and temperature. For large temperatures, QD is strong 
and the degree of mixedness is high, while for zero temperature the asymptotic final 
state is pure. With increasing squeezing parameter and initial correlation, QD becomes 
stronger, but the asymptotic value of the degree of QD does not depend on the initial 
squeezing and correlation, it depends on temperature only. 

(2) Using the criterion of CC, we have shown that the general tendency is that 
CC are less and less strong with increasing time and temperature. For a long enough 
interval of time we have a significant degree of CC, but at t — > oo there are no CC in the 
case of an asymptotic Gibbs state. With increasing squeezing parameter, CC become 
stronger, but the variation of the correlation coefficient has a small influence on the 
behaviour of the degree of CC. 

(3) During a finite interval of time the system interacting with the thermal bath 
manifests simultaneously both QD and CC, so that a true quantum to classical transi- 
tion takes place and the system recovers classicality in a significant measure. CC are 
expressed by the fact that the Wigner function has a peak which follows (exactly for 
A = fi and approximately for A ^ fi) the classical trajectory in phase space and QD 
is expressed by the loss of quantum coherence in the case of a thermal bath at finite 
temperature. For an initial Gaussian quantum state, Wigner function is positive for all 
times, so that it represents a true classical probability distribution in phase space. 

(4) The expressions of the degree of QD and CC confirm the relative competition 
which appears between QD and existence of CC, since decoherence implies a spreading of 
the Wigner distribution function along the trajectory in phase space, whereas CC require 
the existence of sharp peaks in the Wigner function. However, there exists a broad 
compromise regime in which QD and CC can hold well simultaneously. Consequently, 
classicality is a temporary phenomenon, which takes place only at some stages of the 
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dynamical evolution of the system, during a definite interval of time. 

(5) We determined the general expression of the decoherence time, which shows 
that it is decreasing with increasing dissipation, temperature and squeezing. We have 
also shown that the decoherence time has the same scale as the time after which sta- 
tistical fluctuations become comparable with quantum fluctuations, as expected, and 
the values of these scales become closer with increasing temperature and squeezing. Af- 
ter the decoherence time, the decohered system is not necessarily in a classical regime. 
There exists a quantum statistical regime in between. For the considered open system at 
a finite temperature, the uncertainty relation ( 15^1) holds, which interpolates between the 
Heisenberg relation at zero temperature (1551) and the high temperature classical statis- 
tical relation fl56|) . Only at a sufficiently high temperature, when the spin statistics can 
be represented by the Maxwell-Boltzmann distribution, the system can be considered 
in a classical regime. 

The study of classicality using QD and CC leads to a deeper understanding of the 
quantum origins of the classical world. As a result of the progress made in the last 
two decades, the quantum to classical transition has become a subject of experimental 
investigations, while previously it was mostly a domain of theory [2j [3] . The issue of 
quantum to classical transition points to the necessity of a better understanding of 
open quantum systems. The Lindblad theory provides a self-consistent treatment of 
damping as a general extension of quantum mechanics to open systems and gives the 
possibility to extend the model of quantum Brownian motion. The obtained results 
in the framework of the Lindblad theory are a useful basis for the description of the 
connection between uncertainty, decoherence and correlations (entanglement) of open 
quantum systems with their environment. 
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